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ABSTRACT 

We measure apparent velocities (w ap p) of the Ha and H/3 Balmer line cores for 449 non-binary 
thin disk nor mal DA white dwarfs (W Ds) using optical spectra taken for the ESO SN la Progenitor 
surveY (SPY; Nap iwotzki et al.l 120011) . Assuming these WDs are nearby and co- moving, we correct 
our velocities to the Local Standard of Rest so that the remaining stellar motions are random. By 
averaging over the sample, we are left with the mean gravitational redshift, (v g ): we find (i> g ) = 
( u ap P ) = 32. 57± 1.17km s _1 . Using the mass-radius relation from evolutionary models, this translates 
to a mean mass of 0.6471q'q}4 M Q . We interpret this as the mean mass for all DAs. Our results are 
in agreement with previous gravitational redshift studies bu t are significantly higher than all previous 
spectroscopic determinations except the recent findings of iTremblav fc Bergeron! (|2009ft . Since the 
gravitational redshift method is independent of surface gravity from atmosphere models, we investigate 
the mean mass of DAs with spectroscopic T e g both above and below 12000 K; fits to line profiles give 
a rapid increase in the mean mass with decreasing T e g. Our results are consistent with no significant 
change in mean mass: (M) hot = 0.640 ± 0.014 M and (M) c ° o1 - n Rs«+ 035 i 
Subject headings: 



0.686™ M . 



1. INTRODUCTION 

Nearly all stars end their trek through stellar evolution 
by becoming white dwarfs (WDs). Hence, properties of 
WDs can provide important information on the chemi- 
cal and formation history of stars in our Galaxy, as well 
as on the late stages of stellar evolution. Of these stel- 
lar properties mass is one of the most fundamental, and 
though there are several methods for mass determination 
of WDs, each has its limitations. 

The most- widely used WD mass determination method 
involves comparing predictions from atmosphere models 
with observations to obtain effective temperatures (T e g) 
and/or surface gravities (\ogg). One can then compare 
these qu antities with predictions from evolutionary mod- 
els (e .g., lAlthaus fc Benvenutolll998t Montgomery et al.l 
1999). Shipmanl (| 1979( 1. IKoester et al.l l| 19791 ). and 
McMahan (1989|) use radii determined from trigonomet- 
ric parallax measurements along with T e g from pho- 
tometry to determine masses. Of course this tech- 
nique is limited to target stars with measured paral- 
laxes, so users of photometry have more often used ob- 
served color indices to determine both T R ff and log g (e.g. , 
Koester et aDl979UWegnen|1979HShipman fe Sasslll98Clt 
Weidemann fe Koesterlll984l : lFontaine et al.ll985l). Wit h 
the exception of the parallax variant (jKilic et al.l 12008). 
the photometric method is seldom used in recent WD 
research. 

Another variant of this method uses mainly spec- 
troscopic rather tha n photometric obse r vations (e.g., 
Bergeron et"aT1 Il992t iFinlev et all 119971 : iLiebert et all 
20051) . With more recent large-scale surveys, such as 
SPY (s ee Section [5]) and the Sloan Digital Sky Survey 
(SDSS; I York et alT^OOOl) . the comparison of observed 
WD spectra with spectral energy distributions of theo- 
retical atmosphere models has become the primary WD 
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mass determination m ethod, yielding masses for a large 
number of WDs (e.g., IKoester et all 120011 : iMadei et all 
2004| IKeoler et al.ll2007D . 

When applied to cool WDs (T eS < 12000 K), however, 
the reliability of this primary method breaks down: a sys- 
tematic increase in the mean log g for DAs with lower T P ff 
has repeatedly shown up in analyses (e.g. . ILiebert et all 
120051 iKepler et all 120071 : iDeGennaro et all 120081). This 
"log g upturn" , discussed thoroughly by Bergero n et all 
(2007) and by IKoester et all (|2009aD . is generally be- 
lieved to reflect shortcomings of the atmosphere models 
— specifically our understanding of the line profiles — 
rather than a real increase in mean mass with decreasing 
T e ff as an increasing mean log g would imply. 

Other mass determination methods that are indepen- 
dent of atmos phere models include the astr ometric tech- 
nique (e.g., iGatewood fc Gatewoodl I1978D an d pulsa- 
tional mode analysis (e.g.. iWinget et al.lll991ll . Unfor- 
tunately, neither of these methods are widely applicable 
to WDs. The former requires stellar systems with multi- 
ple stars, and the latter is limited to WDs and pre-white 
dwarfs which lie in narrow T e ff ranges of pulsational in- 
stability. 

Another method that is mostly atmosphere model- 
independent uses the gravitational redshift of absorp- 
tion lines; this is the one that will be the focus of this 
paper. The difficulty in disentangling the stellar radial 
velocity shift from the gravitational redshift has caused 
this method to only be used for WDs in common proper 
motion binaries or o pen clusters (iGreenstein fc Trimble! 
19671: IKoester! fl987l : IWegner fc Reidl 119911 : IReidl 119961 : 
Silvestri et al.l 120011 ). The simplicity of this method, 
however, prompts us to extend the investigation beyond 
those cases. 

In this paper, we will make two main points: (1) by 
using a large, high-resolution spectroscopic dataset, we 
can circumvent the radial velocity-gravitational redshift 
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Fig. 1. — Distribution of targets in Galactic longitude I and latitude b. We mark the targets in our main sample as black points 
and the thick disk WDs as pink squares. We indicate the direction of the movement of the Sun with respect to the LSR (blue cross; 
IKerr fc Lvn dcn-Bcll 1986). Since the observations are from the ESO VLT in the Southern Hemisphere, no targets with a declination above 
+30° are in our sample, hence the gap in the left side of the plot. 



degeneracy to measure a mean gravitational redshift of 
WDs in our sample and use that to arrive at a mean mass; 
and (2) since the gravitational redshift method has the 
advantage of being independent of surface gravity from 
atmosphere models, we can use it to reliably probe cool 
DAs (T e ff < 12000 K), thus providing important insight 
into the "log g upturn problem" as groups continue to 
impr ove upon those models (e.g.. iTremblav fc Bergeron! 
[2001 . 

2. GRAVITATIONAL REDSHIFT 

In the weak-field limit, the general relativistic effect of 
gravitational redshift (z) can be understood, classically, 
as the energy (E) lost by a photon as it escapes a gravi- 
tational potential ($) well: 



-AE 
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il) 



The fractional change in energy can be rewritten as a 
fractional change in observed wavelength (-AE/E = 
AA/A). In our case, the gravitational potential is at the 
surface of a WD of mass M and radius R. In terms of a 
velocity, the gravitational redshift is 
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where G is the gravitational constant, and c is the speed 
of light. 

For WDs, Vg is comparable in magnitude to the stellar 
radial velocity v T , both of which sum to give the apparent 
velocity we measure from absorption lines: t> ap p = v g +v r . 
These two components cannot be explicitly separated for 



individual WDs without an independent v T measurement 
or mass determination. 

The method of this paper is to break this degener- 
acy not for individual targets but for the sample as a 
whole. We make the assumption that our WDs are a 
co-moving, local sample. After we correct each w app to 
the Local Standard of Rest (LSR), only random stellar 
motions dominate the dynamics of our sample. We as- 
sume, for the purposes of this investigation, that these 
average out. Thus the mean apparent velocity equals the 

) = ( v g)- The idea of 
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mean gravitational redshift: 

averaging over a group of WDs to extract a mean gravita - 
tional redshift is not new (|Greenstein fe Trimble! [1967). 
but the availability of an excellent dataset prompted its 
exploitation. We address the validity of the co-moving 
approximation in Section 14.11 



3. OBSERVATIONS 

We use spectroscopic data from the European South- 
ern Observatory (ESQ) SN la Progenitor surveY (SPY; 
Napiwotz ki et al.l [2001). These observations, taken 
using the UV- Vi sual Echelle Spectrograph (UVES; 
iDekker et al.l [20001 at Kueyen, Unit Telescope 2 of the 
ESO VLT array, constitute the largest, homogeneous, 
high-resolution (0.36 A or ~ 16 km s _1 at Ha) spectro- 
scopic dataset for WDs. We obtain the pipeline-reduced 
data online through the publicly available ESO Science 
Archive Facility. 

3.1. Sample 

As explained in iNapiwotzki et all ()200lD , targets for 
the SPY sample come from: the white dwarf catalog of 
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Wisotzki et all 120001 iChristlieb et~afl 1200 ID. the Ham- 
burg; Quasar Survey ( Hagen et al.l Il995t iHomeier et al.l 



1998) , the Montreal- C ambridge- Tololo survey (MCT; 
Lamontagne et all I2000T) . and the Edinburgh-Cape sur- 
vev (EC: Kilkenny et all 119971 ). The magnitude ol the 
targets is limited to B < 16.5. 

Our main sample consists of 449 analyzed hydrogen- 
dominated WDs (see Figure [T] for the distribution of 
targets in Galactic coordinates). This is the subset of 
the SPY sample that meets our sample criteria (ex- 
plained below) and that shows measurable v app in the Ha 
(and H/3) line cores while not showing measurable v app 
variations. A variable velocity across multiple epochs 
of observation suggests binarity. The method of SPY 
to search for double degenerate systems is to detect 
variable radial velocity. For our study, however, we 
are interested only in non-binary WDs since these pre- 
sumably have no radial velocity component in addition 
to random stellar motion after being corrected to the 
LSR. We exclude known double degenerates and common 
proper motion binary systems (iFinlev fc Koesterl 119971: 
Jordan et al.lll998tlMaxted fc Marshlll999):IMaxted et al l 
20001: ISilvestri et all 120011: IKoester etaflfeoogbfl even if 
we do not find them to show variable v app . 

We choose "norm al" DAs (criterion 1) from 
IKoester et al.l (|2009bl ). Classification as a normal 
DA does not include WDs that exhibit He absorption 
in their spectra in addition to H absorption, and it does 
not include magnetic WDs. In a subsequent paper we 
will investigate the sample of 20 helium-dominated WDs 
for which we observe H absorption. 

For our main sample, we are also only interested in 
thin disk WDs (criterion 2), so we exclude ha lo and thick 
disk candidates as kinema tically classified bv lPauli et al.l 
(poM) and lRichter et~aT1 $mm\ . We assume the rest are 
thin disk objects, the most numerous Galactic compo- 
nent. Our sample selection is also con sistent with the re - 
sults for the targets in common with iSion et al.l (|2009l) . 
iRichter et all (|2007l ) find only 2% and 6% of their 632 
DA WDs from SPY to be from the h alo and thick disk , 
respectively. For WDs within 20 pc, ISion et aTl (|2009f ) 
find no evidence for halo objects and virtually no thick 
disk objects. We note that unique identification of pop- 
ulation membership for WDs is difficult and often not 
possible because of ambiguous kincmatical properties. 
B ased on corr e ctions for these intrinsic contaminations 
bv lNapiwotzkil <|2009f ) , we expect any residual contamina- 
tion in our sample to be at most ~ 6%. A contamination 
this size will have a negligible impact on our conclusions. 
We explain the significance of requiring thin disk WDs in 
Section 14. 1| and we explore a mini-sample of thick disk 
WDs in Section I5T41 

The gravitational redshift method becomes very diffi- 
cult for hot DAs with 50000 K > T cff > 40000 K (see the 
T e ff gap in Figure \5§ . As the WD cools through this T e g 
range, the Balmer line core, which we use to measure v app 
(Section |4j, disappears as it transitions from emission to 
absorption; fortunately only ~ 5% of the DAs from SPY 
lie in this range. 

4. VELOCITY MEASUREMENTS 

In the wings of absorption lines, and in particular, 
for the hydrogen Balmer series, the effects of colli- 
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Fig. 2. — Example UVES spectrum for a target in our sample. 
We measure v ap p by fitting Gaussian profiles (solid, red lines) to 
the non-LTE Balmer line cores using a non-linear least-squares 
fitting routine. The line cores are well-resolved, allowing for precise 
centroid determinations. 

sional broadening cause asym metry, making it difficult 
to measure a velocity c entroid (jShipman fc Mehanll 19761 : 
iGrabowski et al.lll987l) . These effects are much less sig- 
nificant, however, in the sharp, non-LTE line cores, and 
furthermore with decreasing principal quantum number, 
making both the Ha and H/3 line cores suitable options 
for measuring an apparent velocity i> app - Higher order 
Balmer lines are intrinsically weaker (the H7 line core, 
for example, is seldom observable in our data), so finite 
S/N prevents the number of observable H/3 line cores 
from matching the number of observable Ha line cores. 

We measure v app for each target in our sample by fit- 
ting a Gaussian profile to the Ha line core using GAUSS- 
FIT, a non-linear least-squares fitting routine in IDL (see 
Figured] for an example). When available, we combine 
this measurement with that of the H/3 line core cen- 
troid as a mean weighted according to the uncertain- 
ties returned by the fitting routine. We include H/3 line 
core centroid information in 372 of our 449 u app mea- 
surements. If multiple epochs of observations exist, we 
combine these measurements as a weighted mean as well. 
Apparent velocity measurements of a given observation 
(i.e., Ha and H/3 line core centroids) are combined before 
multiple epochs. 

Table [1] (full version available on-line) shows our mea- 
sured v app for Ha and H/3 (when observed) for each ob- 
servation. 

4.1. Co- Moving Approximation 

We measure a mean gravitational redshift by assuming 
that our WDs are a co-moving, local sample. With this 
assumption, only random stellar motions dominate the 
dynamics of our targets; this falls out when we average 
over the sample. 

For this assumption to be valid, at least as an approx- 
imation, our WDs must belong to the same kinematic 
population; in the case of this work, this is the thin disk. 
We achieve a co-moving group by correcting each mea- 
sured v app to t he kinematical LSR describ ed by Standard 
Solar Motion (jKerr fc Lvnden-Belllll986l ). 

There are reasons to believe that the targets in our 
sample will not significantly lag behind our choice of 
LSR due to asymmetric drift. Although WDs are con- 
sidered "old" since they are evolved stars, it is the to- 
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Fig. 3. — Distri bution of d istances (from spectroscopic parallax) 
of SPY WDs from Pauli et al. ( 2006) . The shaded, green histogram 
shows the targets in our sample. The mean is 94.5 pc; the median is 
89.2 pc. These distances are short enough to support our co-moving 
approximation. We list the number of targets in each distribution. 

tal age of the star (main sequence lifetime r nuc and 
cooling time r CO oi) that is of consequence. WDs with 
M ~ 0.6 M ^ have main sequence progenitors with M ~ 
2 Mq (e.g.. iWilliams et all 120091) This corresponds to 
Tnuc ~ 1.4 Gyr (|Girardi et al . 1 120001) . r CO oi is on the order 
of a few hundred million years for most of the WDs in 
our sample (T e g of a few times 10 4 K) and ~ 2.5 x 10 9 yr 
for our coolest WDs (T c g ~ 7000 K); the total age spans 
a range of roughly 1.5 to 4 Gyr (F/G type stars). 

We also make certain that our WDs reside at distances 
that are small when compared to the size of the Galaxy, 
thereby making systematics introduced by the Galactic 
kinematic structure negligible. Figure [3] shows the d is- 
tances (from spectroscopic parallax; iPauli et al1l2006| ) to 
the targets in our sample. The mean distance of the 
targets in the histogram is less than 100 pc, and all are 
within 600 pc. Over these distances, the velocity disper- 
sion with varying height a bove the disk remains modest 
(|Kuiiken fc Gilm ore 198jJ), and differential Galactic ro- 
tation changes very little as well (~ 3 km s^ lFich et al.l 
1989). In Section [5.3.21 we perform an empirical check 
to the assumptions made in this section. 

5. RESULTS 

5.1. Mean Apparent Velocities 

We present the distribution of our measured apparent 
velocities in Figure [U Table [T] lists individual apparent 
velocity measurements, and mean apparent velocities are 
in Table H 

Though our main method uses information from both 
the Ha (Column 7 of Table [TJ and H/3 (Column 9) line 
cores to determine t; ap p for a given observation (Column 
11), we also perform our analysis using Ha only and H/3 
only. We measure H/3 line core centroids for 382 of our 
449 targets. 

Figure [5] shows the distribution of measurement un- 
certainties associated with each target. H/3 centroid de- 
terminations are typically less precise than those for Ha 
(see Column 6 of Table [3]) , which is expected since the 
Ha line core is nearly always better-defined. We find 
that the improved precision achieved by combining Ha 
and H/3 information is not significant when determining 
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Apparent Velocity (km s ) 

Fig. 4. — Histograms of measured apparent velocities ri a pp with 
a bin size of 5 km s" 1 . The mean t> app for all targets in our sample 
(shaded) is 32.57 ± 1.17km s — 1 ; the median is 31.94km s~ ; the 
standard deviation is 24.84 km s _1 . Using t> app measured from Ha 
only (red, descending lines): the mean n app is 32.69 ± 1.18km s _1 ; 
the median is 32.05 km s _1 ; the standard deviation is 24.87 km 
s _1 . Using v app measured from H/3 only (blue, ascending lines): 
the mean n app is 31.47 ± 1.32 km s — 1 ; the median is 31.55 km s _1 ; 
the standard deviation is 25.52 km s _1 . The overplotted curves are 
the Gaussian distribution functions used to determine Monte Carlo 
uncertainties. We list the number of targets in each distribution. 

the uncertainties to our mean apparent velocities. These 
uncertainties are dominated by sample size. In fact, we 
must increase (worsen) our typical measurement error of 
~ 2 km s _1 to ~ 10 km s _1 to notice a ~ 7% increase 
in the error of the mean; a monstrous leap to measure- 
ment uncertainties of ~ 50 km s _1 enlarges the error of 
the mean by a little more than a factor of 2. Thus, using 
Ha (or H/3) centroids only is sufficient for the kind of in- 
vestigation employed in this paper, and lower resolution 
observations are also suitable as long as the Balmer line 
core is resolved. 

The quoted uncertainties of the mean apparent veloci- 
ties (Column 4 of Table [3]) come from Monte Carlo sim- 
ulations. For each sample, we recreate a large number of 
instances (10000) of the v app distribution by randomly 
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Fig. 5. — Histograms of apparent velocity measurement uncer- 
tainties Svupp corresponding to the samples in Figure [J] The bin 
size is 0.4km s — 1 . The overplotted curves are the empirical dis- 
tribution functions used to determine Monte Carlo uncertainties. 
Notice that measurements of the Ha line core are more precise 
than for H/3. For aesthetics, we leave off two H/3 <5i; app of 13.06 
and 17.57 km s _1 . 
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Fig. 6. — Dist ribution of spectrosco pically determined T e ff of 
normal DAs from Kocstcr ct al. (2009b) (dashed, black histogram). 
The bin size is 1500 K. The solid, green histogram shows the non- 
binary thin disk SPY targets, and the shaded, purple histogram 
shows the targets in our sample. The mean is 19400 ± 300 K; the 
median is 17611 K.; the standard deviation is 9950 K. The ovcrplot- 
ted curve is the empirical distribution function used to determine 
Monte Carlo uncertainties. We list the number of targets in each 
distribution. 

sampling from a convolution of the empirical i> app dis- 
tribution (Gaussian characterized by the parameters in 
Columns 3 and 5 of Table [3]) and the empirical measure- 
ment error distribution. We adopt the standard devia- 
tion of the resulting simulated mean values as our formal 
uncertainties. Since the input distributions for our sim- 
ulations are empirical, our uncertainties are subject to 
the normal limitations of Frequentist statistics. We plot 
the empirical distribution of our main sample in Figure 
[4] (black curve) along with the distributions for the Ha 
(dashed, blue curve) and H/3 (dashed, red curve) sam- 
ples. The corresponding empirical distributions of our 
measurement uncertainties are in Figure [5] 

For convenience, Table[3]also lists the quantity (M/R), 
which is propor tion al to (v g ) (equation [5]) and, as we 
argue in Section |4~T1 (v a pp)- 

5.2. Mean Masses 

The mean apparent velocity {v app ) (or (M/R)) is our 
fundamental result since it is this quantity that is model- 
independent. To translate this to a mean mass (Table 
HJ, we must invoke two dependencies: (1) we need an 
evolutionary model to give us a mass-radius relation, and 
(2) since the WD radius does slightly contract during its 
cooling sequence, we need an estimate of the position 
along this track for the average WD in our sample (i.e., 
a mean T g). 

Our evolutionary models use M-Rc/M* = 10 -2 and 
Mh/M+ = 10~ 4 for the surface-layer masses; these are 
canonical values deriv e d from evoluti onary studies (e.g., 
lLawlor fe MacDonaldl [2006). See Montgomery et al.l 
(|1999| ) for a more complete description of our models. 
Our dependency on evolutionary models is small. We are 
interested in the mass-radius relation from these mod- 
els, and this is relatively straight-forward since WDs are 
mainly supported by electron degeneracy pressure, mak- 
ing the WD radius a weak function of temperature. We 
estimate that varying the C/O ratio in the core affects 
the radius by less than 0.5%, whereas changing M-r/M* 
from 10 -4 to 10 -8 results in about a 4% decrease in 



radius. See Section 15.3.11 for more discussion on the de- 
pendency of the hydrogen layer mass. 

Figure[7]plots M/R versus T e g with cooling tracks from 
evolutionary models for a range of WD masses. We use 
(T off ) = 19400 ± 300 K from the spect roscopically de- 
termined values of Kocstcr ct al. (2009b) (see Figure [6]) , 
and, after plotting (M/R) from Table 02 we interpolate 

to arrive at a mean mass of 0.647^q'q}4 M© for 449 non- 
binary thin disk normal DA WDs from the SPY sample. 

To compare this result with that of the spectroscopic 
meth od, we use atmospher ic parameters logg and T e g- 
from iKoester et alJ (2009b) along with the mass-radius 
relation from evolutionary models to der ive individual 
masses for 441 of the targets in our sample (Kocstcr ct al. 
(2009b]) do not publish individual WD masses). We de- 
rive a sharply peaked mass distribution (Figure [8]) with 
width (not uncertainty) a = 0.13 M Q and a mean mass 
of 0.575 ± 0.002 M — significantly lower than the value 
we obtain from the gravitational redshift method. We 
compute the error of the mean using Monte Carlo sim- 
ulations following the same method described in Section 
I5.1l except instead of using a single Gaussian to represent 
the mass distribution, we use multiple Gaussians (curve 
in Figured]). 

5.3. Systematic Effects 
5.3.1. From Evolutionary Models 

The hydrogen layer mass in DAs is believed to be in 
the range of 10" 4 > M H /M* > 10" 8 , constrained by hy- 
drogen shell burning i n the late stages of stellar evo lution 
([Althaus et al.l 120021: lLawlor fc MacDonaldl I2006D and 
convective mixing ((Fontaine fe Wesemaellll997fl. I n their 
asteroseismological studies. iBis choff-Ki m et al.1 (2008) 
also find evidence to support this range of hydrogen 
layer masses, and this is co nsistent with the results of 
iCastanheira fc Keplerl (|2009ft . 

Our evolutionary models use the fiducial value of 
M-r/M+ = 10 ~ 4 for "thick" hydrogen layers. First, 
this is suggeste d by the pre- white dwarf ev olutionary 
models of, e.g., lLawlor fc MacDonal d (2006), who find 
that the overwhelming majority of their DA models 
have thick hydrogen layers. Second, if thin layers were 
the norm, then convective mixing below 10000 K would 
lead to a disappearance of DAs at these temperatures 
(Fontaine & Wesemacl 1997). Both of these reasons lead 
us to choose thick hydrogen layers for our models. 

We find that using a midrange hydrogen layer mass 
of Mn/M+ = 10 -6 decreases the mean mass we derive 
for our main sample by 0.012 M Q , while using a thin 
layer mass of M^/M^ = 10~ 8 decreases the derived mean 
mass by an additional 0.003 M Q (total mass difference of 
0.015M Q ). Assuming no hydrogen layer (Mh/M* = 0) 
yields a mean mass that is ~ 0.018 M Q lower than that 
obtained with the fiducial value of Mjj/A/* = 10~ 4 . 

It is worth noting that the spectroscopic method shares 
this dependency on evolutionary models and that most 
of the studies listed in Table El i ncluding iLiebert et al.1 
2005), iKepler et all ([20071 ) and iTremblav fc Bergeron! 



2009), employ mass-radius relations that use thick hy- 



drogen layers. Column 7 of Table [5] notes the assumed 
hydrogen layer mass in the evolutionary models used in 
each study. Furthermore, our results are qualitatively 
less sensitive to the mass-radius relation: for the grav- 
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Mass/Radius (M /R o ) 

Fig. 7. — {left) Plot of M/R versus T a g with cooling tracks from evolutionary models for a range of WD masses. The intersection of the 
mean measured apparent velocity u a pp {vertical, black line) and mean T e fj from Figure [6] {horizontal, purple line) indicates a mean mass 
of 0.647^q'q^ Mq. {right) A version of Figure [5] with an abbreviated temperature range. We leave off 13 WDs with T e g > 50000 K from 
the plot. 
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Fig. 8. — Distribution of spectroscopic masses of normal DAs from Kocstcr et al. (2009b) we derive using the published atmospheric 
parameters logg and T a g {dashed, black histogram). The bin size is 0.025 Mq. The solid, green histogram shows the non-binary thin disk 
SPY targets, and the shaded, orange histogram shows the targets in our sample. The means are 0.567 ± 0.002 Mq {vertical, black line), 
0.580 ± 0.002 Mq {vertical, green line), and 0.575 ± 0.002 Mq {vertical, orange line), respectively. Notice that the mean spectroscopic 
masses are similar, indicating that the application of our sample criteria to SPY is not introducing additional systematic effects. All the 
means are also significantly less than the mean mass derived from the gravitational redshift method {vertical, blue line). The overplotted 
curve is the empirical distribution function used to determine Monte Carlo uncertainties. We list the number of targets in each distribution. 
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itational redshift method, v g oc M/R, while the sur- 
face gravity used by the spectroscopic method scales as 
g oc M/R 2 . 

5.3.2. Dynamical 

We use the kinematical LSR describ ed by Standard 
Solar Motion (jKerr fc Lvnde n-Bcll 19861) as our reference 
frame for the co-moving approximation. To determine if 
this is a suitable choice, we investigate (i> app ) in the U , V 
or W directions (by convention, U is positive toward the 
Galactic center, V is positive in the direction of Galactic 
rotation, and W is positive toward the North Galactic 
Pole). 

For 237 targets in the direction of the Galactic center 
(I < 90° or I > 270°) and 212 opposite the Galactic 
center (90° < I < 270°), (v app ) = 31.81±1.71 and33.43± 
1.64 km s~ x , respectively. In the direction of the LSR 
flow (/ = 90°, 6 = 0°; 196 targets) and opposite the flow 
(253 targets), (v app ) = 33.61 ± 2.09 and 31.77 ± 1.34km 
s -1 . North (185) and south (264) of the Galactic equator, 
( Uapp ) = 31.59 ± 1.84 and 33.26 ± 1.53 km s" 1 . 

These empirical checks provide independent evidence 
that the local WDs in our sample move with re- 
spect to kinematical LSR with the following values: 
(U, V, W)=(-1.62 ± 3.35, +1.84 ± 3.43, -1.67 ± 3.37) km 
s _1 , which is consistent with no movement relative to the 
LSR. Therefore, we find our choice of reference frame to 
be suitable for this study. 

5.3.3. Observational 

SPY targets are magnitude-limited to B < 16.5, but 
these targets come from multiple surveys with varying 
selection criteria, maki ng the combined crite ria difficult 
to precisely determine (|Koester et al.l [200 9b ). For this 
reason, our results pertain mostly to non-binary thin disk 
normal DA WDs from SPY. Although the selection bias 
is likely to have a minimal effect, a detailed comparison 
of our results with that of the general DA population 
awaits a closer examination o f the selection criteria (see 
iNapiwotzki et~aI|[200lT [200l . 

If we approximate our sample to be free of any target 
selection bias, our crude estimates show that we have a 
net observational bias toward lower mass WDs. There 
are two competing effects: first, at a given T e g, a larger 
mass (smaller radius) results in a fainter WD, thus bias- 
ing the detection of fewer higher mass WDs over a given 
volume, and second, a larger mass (smaller radius) also 
results in a slower cooling rate due to a larger heat ca- 
pacity as well as a diminished surface area. This means 
more higher mass WDs as a function of T e g . We estimate 
the observational mass bias correction as follows: 

Let P(M) be the distribution of WDs as a function 
of mass for a magnitude-limited sample of WDs. For 
simplicity, we take it to have the form of a Gaussian; we 
take the mean to be (M) ~ 0.65 M© and a ~ 0.1 M . 
As a reference, the spectroscopic mass distribution of 
DAs shows a sharp Gaussian-like peak with high and 
low mass wings (e.g.. iBergeron et alj|1992t fLiebert et al.l 
[20051 IKepler et al.ll2007D 

Effect (1): Ignoring color, the apparent flux of a star 
scales as -F app ~ L+/D 2 and the luminosity as ~ 
i? 2 T c 4 ff , where L*, R and T e g are the luminosity, radius, 
and effective temperature of the star; D is its distance. 



In the non-relativistic limit, the radius R of a WD scales 
as R oc Al" 1 / 3 (jChandrasekharl 1 1 9390 . and for a (mod- 
erately relativistic) 0.6 M© WD this relation is approxi- 
mately R oc M -1 / 2 , so 



M 



(3) 



If -F cu toff is the lower limit on flux for the survey, a given 
WD is visible out to a distance of 



D 



F, 



cutoff 



1/2 



cx 



MV2- 



(4) 



If we make the simplifying assumption that all the WDs 
are at the observed average temperature (T c g) and that 
they are distributed uniformly, the volume V in which a 
WD is visible is 



V ~ D 3 oc M~ 



(5) 



Thus, P(M) is biased by this factor. 

Effect (2): From simple Mestel theory (|Mestellll952l) . 
the WD cooling time r scales as 



r oc 



M\ 7 
17 



which, from equation |3l yields 



1 eS 



(6) 



(7) 



Again, assuming that the WDs are all at (T c s), the ob- 
served distribution will be biased by a factor of r oc 
M w ' 7 . 

Thus, the final biased distribution we observe is given 
by the product of these factors: 



P hias {M)^VTP{M) 

(xAr 1/u P(M). 



(8) 



This very weak mass bias results in {M)bj a s = 0.649 M©, 
which is a mass bias of AM — —0.001 M©. While this is 
just a crude estimate, it suggests that the bias correction 
is likely much smaller than the size of our stated random 
uncertainties. 

5.3.4. Mass Conversion 

In our mean mass determination in Section [5. 2 [ we im- 
plicitly assume that (M/R) = (M)/(R). These quanti- 
ties are not entirely equal, and by performing an estimate 
using a simple analytical form for the WD mass distri- 
bution, we find that there is a difference of ~ 0.5% (i.e., 
{M/R) ~ 1.005 x (M)/{R)), which we consider to be a 
negligible systematic. 

5.4. Thick Disk DAs 

The kinematics of thick disk stars prohibit us from 
placing them in the same co-moving reference frame as 
thin disk stars. In Section [5.3. 2[ we show that the kine- 
matical LSR described by Standard Solar Motion is a 
suitable choice of reference frame for the SPY thin disk 
WDs. As expected, using w app of our thick disk targets 
corrected to that LSR (the reference frame suitable for 
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Fig. 9. — Distribution of T e g versus logg for 41 9 of our WDs. 
Spectro scopic parameters for all targets are from IKoester et al.l 
(2009b). Notice the abrupt increase in the mean log q around 
12000 K. We also plot cooling tracks from evolutionary models for 
0.5, 0.6 and 0.7 Mq WDs. 

the thin disk) give discrepant values for (i> app ) in oppo- 
site directions. Since our thick disk sample is small (26 
targets), our (v app ) uncertainties are too large to discern 
a suitable reference frame. If we correct by the average 
lag in rotational velocity of t he thick disk with r espect 
to the thin disk (~ 40 km s" 1 ; iGilmore et al.lll989l ). then 
(^app) = 32.90 ± 9.59 km s _1 for our thick disk sample. 
Individual i; app measurements are listed in Table [2] Us- 
ing (T cff ) = 19960K, we find (M) = 0.652±g;?^ M Q , 
which is evidence that the mean mass of thick disk DAs 
is the same as for thin disk DAs. 

One should also notice that the dispersion of v app (Col- 
umn 5 of Table [3]) is clearly larger than that for the thin 
disk DAs. Since the v app distribution is a convolution 
of the true mass distribution and the random stellar ve- 
locity distribution, this is consistent with a larger veloc- 
ity dispersion as exp ected for the thick disk population 
(IGilmore et al.l[l989h . 

6. DISCUSSION 
6.1. The Logg Upturn 
6.1.1. The Problem 

A major problem plaguing the field of WDs is the 
apparent systematic increase in mean log g for DAs to- 
ward low (< 12, 000 K) T e ff, as determine d from spectro- 
scopi c fitting of absorption line profiles (jBergeron et al.l 
120071 IKoester et~aLl |2009aD . This increase is absent 
in photometric logg determin ations (|Kepler et al.l 120071 : 
lEngelbrecht fc Koe stcr 2007), which are not strongly 
dependent on line profiles. A number of effects are 
known to exist that make theoretical line profile model- 
ing for cool WD atmospheres more difficult than for hot- 
ter WDs, such as helium contamination from dred ge-up 
(jBergeron et al.l 119901: iTremblav fc Bergeron! I2008D and 
the tr eatment of convective efficiency (jBergeron et al.l 
Il995aj ). Neither of th ese, however, seem t o solve the 
log g upturn problem (|Koester et al.l l2009af l . and since 
no strong hypotheses have been put forth to explain a 
real increase in mean mass (|Kepler et al.ll200"7T) , the fault 
most likely lies with the atmosphere models or with the 
limitations of these models. 

The number of studied cool WDs is already rela- 



FlG. 10. — Histogram of measured apparent velo cities u a pp for 
targets with spectroscopically determined T e g from IKoester et al.l 
(2009b). The bin size is 7.5 km s" 1 . The green histogram with 
ascending lines corresponds to targets with T a g > 12000 K and the 
pink histogram with descending lines to 12000 K > T c ff > 7000 K. 
The mean i> app for the green histogram is 31.61 ± 1.22 km s _1 ; 
the median is 31.71km s ; the standard deviation is 23.22 km 
s _1 . The mean t) app for the pink histogram is 37.50 ±3.59 km s -1 ; 
the median is 36.20km s~ x ; the standard deviation is 31.00 km 
s _1 . The overplotted curves are the Gaussian distribution func- 
tions used to determine Monte Carlo uncertainties. We list the 
number of targets in each distribution. 

tively low due to the inherent difficulty of observing 
cool objects, but the addition of the log g upturn prob- 
lem and the subtleties of cool WD atmosphere model- 
ing has thus far kept that number low by prompting 
many spectroscopic analyses to be design e d to exclude 
coole r WDs (e.g., iBergeron et al.l 119921: iMadei et all 
l200l iLiebert et all 1200a Kepler et al.l 120071) . This is 
tremendously unfortunate. Understanding cool WDs has 
broad astrophysical relev ance, such as in dete rmining the 
age of the Galactic disk (jWinget et al.lll987l ) and in set- 
ting constraints on the physics of cry stallization in high- 
density plasmas (jWinget et al.1120091) . 

Furthermore, decades of focus on hotter WDs (due to 
the much larger dataset and due to the neglect of cooler 
WDs) have perhaps given researchers in our field a false 
comfort with these objects. There is a feeling that since 
hot WD atmospheres are more straightforward to model 
than cool atmospheres, the spectroscopic surface gravi- 
ties (and masses) must be correct for the hot WDs and 
not for the cool WDs, given the logg upturn problem. 
Recent improved calculations for S tark broadening of hy- 
droge n lines in WD atmospheres (jTremblav fc Bergeron! 
2009) show that hot WD modeling is still maturing. 

6.1.2. Avoiding the Upturn 

The gravitational rcdshift method is independent of 
log g determinations from atmosphere models and allows 
us to constrain changes in mean masses across T c g bins. 

Figure M plots sp ectroscopically determi ned values of 
logg and T e ff from IKoester et alj (|2009bD for the tar- 
gets in our sample, clearly exposing the upturn. We 
plot evolutionary models for 0.5, 0.6 and 0.7 M Q DA 
WDs to illustrate how a higher surface gravity implies 
a higher mass and to show the expected weak depen- 
dence on T c g. Using the mass-radius relation from evo- 
lutionary models, we derive mean spectroscopic masses 
{M) hot = 0.563 ± 0.002 M for 358 WDs with T eff > 
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Fig. 11. — Histogram of measured apparent velo cities « app for 
targets with spectroscopically determined T e g from IKoester et al.l 
(2009b). The bin size is 7.5km s _1 . The green histogram with 
ascending lines corresponds to targets with T a g > 12000 K and the 
pink histogram with descending lines to 12000 K > T cff > 7000 K. 
The mean i>app for the green histogram is 31.61 ± 1.22 km s ; 
the median is 31.71km s — 1 ; the standard deviation is 23.22 km 
s — . The mean i) a pp for the pink histogram is 37.50±3.59km S — 
the median is 36.20 km s _1 ; the standard deviation is 31.00 km 
s — ■ The overplotted curves are the Gaussian distribution func- 
tions used to determine Monte Carlo uncertainties. We list the 
number of targets in each distribution. 

Fig. 12. — Similar to Figure[3]but corresponding to targets with 
T off > 12000 K {green) and to targets with 12000 K > T cff > 7000 K 
{pink) . 

12000 K and (M) cooi = 0.666 ± 0.005 M Q for 75 WDs 
with 12000 K > T cff > 7000 K; A{M) = 0.103±0.007M Q . 
The mass differenc e is even larger in the SDSS data; 
Kepler et al] (|2007i ) find (M) hot = 0.593 ± 0.016 M Q and 
(M) c ° o1 = 0.789 ± 0.005 M Q (12000 K> T cS > 8500 K); 
A(M) = 0.196 ± 0.021 M . 



In Figure 1 1 1 1 we show our distribution of 



'app 



(dis- 



tribution of uncertainties in Figure 1121) for targets with 
T c ff > 12000 K (green histogram with ascending lines) 
and with 12000 K> T cff > 7000 K (pink histogram with 
descending lines). The corresponding (w a pp) determina- 
tions are 31.61 ± 1.22 and 37.50 ± 3.59 km s _1 , respec- 
tively, which translates to (M) hot = 0.640 ± 0.014 M Q 
and (Af) c ° o1 = O.686l^M (see Figure H3J). This 
is consistent with no change in mean mass across a 
temperature split at T e g = 12000 K in agreem ent with 
the photometric studie s by Kepler et al.l (|2007l ) and by 
lEngelbrecht fe Koesterl (|2007l ). No previous large spec- 
troscopic study has seen consistency in mean mass across 
these temperatures. 

6.2. Comparison With Previous Studies 

Table [5] lists four studies that employ the gravitational 
redshift method to determine masses for common proper 
motion WDs. Because of the small sample sizes (9, 35, 
34, and 41 WDs), the uncertainties of the mean masses 
found by these studies are relatively large — too large 
to discern a difference in mean mass from that of the 
spectroscopic meth od (ISilvestri et al.ll2001[) . Other than 
with the results of IKoester! (11981 7 whose sample con- 
sisted of only 9 DAs, our mean mass agrees with that of 
all these studies, and we improve upon the uncertainties 
(precision) by more than a factor of two. 

The mean mass of 512 SPY non-binary thin disk nor- 
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Fig. 13. — Same as Figure[7]but for targets with T eff > 12000 K 
{green) and with 12000 K> T eff > 7000 K {pink). (M) hot = 



0.640 ± 0.014 Mq and (M) c ° o1 



,+0.035 - 



mal DAs from IKoester et ahl (|2009bl. as we figure from 
their spectroscopically determined values of log g and 
T G ff, is 0.580 ± 0.002 M Q , and if we restrict the compari- 
son to 441 WDs in our sample, (M) = 0.575 ± 0.002 M . 
Both values are significantly smaller than the mean mass 
we derive using the gravitational redshift method. 

Using atmosphere models tha t implement the new 
Stark broadened line profiles from lTremblay fe Bergeron! 
( 2009) and an updated treatment of the microfield distri- 
bution, the SPY sample shows an increase of ~ 0.03 M 
in the mean mass (Koester, private communication), 
but this resulting mean mass is still significantly less 
than our value. In fact, our mean mass is significantly 
larger than the determinations from all the previous 
spectroscopic studies listed in Table [5] except that of 
ITremblav fe Berger on (2009). 

The recent work of iTremblav fe Bergeron! (|2009f ) uses 
atmosphere models with improved Star k broadening 
calculations to reanalyze the WDs from iLiebert et al.1 
(2005). They find a larger mean mass (0.649 Mq) than 
previously determined for the Palomar-Green sample 
(0.603 M ). The mean mass we derive using the gravita- 
tional redshift method agrees well, thus providing inde- 
pendent observational evidence in support of these im- 
proved atmosphere models. 

7. CONCLUSIONS 

We show that the gravitational redshift method can 
be used to determine a mean mass of a sample of WDs 
whose dynamics are dominated by random stellar mo- 
tions. For 449 non-binary thin disk normal DA WDs 
from SPY, we find (v s ) = (v app ) = 32.57 ± 1.16 km s -1 . 
Using the mass-radius relation from evolutionary models, 
(M) = 0.647l°;o 1 4 M . This is in agreement with the re- 
sults of previous gravitational redshift studies, but it is 
significantly higher than all previous spectroscopic det er- 
minations except that of ITremblav fe Bergeron! p009). 

We find that the targets in our sample move with re- 
spect to the kinematical LSR describe d by Standard So- 
lar Motion (|Kerr fe Lvnden-B cll 1986) with the following 
values: (U, V, W)=(-l.62 ± 3.35, +1.84 ± 3.43, -1.67 ± 
3.37) km s . This is consistent with no movement rela- 
tive to this LSR. 

Our results provide evidence that the mean mass of 



10 



Falcon, Winget, Montgomery, & Williams 



thick disk DAs is the same as for thin disk DAs. 

The gravitational redshift method is independent of 
spectroscopically determined surface gravity from at- 
mosphere models and is insensitive to the logg up- 
turn problem (Section 16. ip . (v app ) — 31.61 ± 1.22 and 
37.50 ± 3.59 km s" 1 for targets with T cff > 12000 K and 
with 12000 K> r c ff > 7000 K, respectively. This trans- 
lates to (Af) hot = 0.640 ± 0.014 M© and (M) c ° o1 = 
0.686lQ 03g M , which disagrees with spectroscopic re- 
sults by showing no significant change in the mean mass 
of DAs across a temperature split at T R ff = 120 00 K. 
This confirms the resuls of iKepler et al.l (|2007l ) and 
lEngelbrecht fc Koesterl (|2007l ) , who find no log g increase 
in their photometric investigations. We are currently ob- 
taining more observations of cool WDs to increase our 
sample size and hence precision of our mean mass deter- 
minations. 
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TABLE 1 

Apparent Velocity Measurements for Normal DA WDs (Abbreviated) 



to 



Adopted LSR Ho H/3 Observation 





^app 




Date 


Time 


Correction 




dv app 


f app 


di? a pp 




app 




Target 


(km s _1 ) 


(km s _1 ) 


(UT1 


(UT1 


(km s _1 ) 


(km. s — ) 


(km s~ ) 


(km s _1 ) 


(km s _1 ) 


(km s — ) 


(km s — ) 




WD 0000-186 


24.530 


0.015 


2000.09.16 


04:53:07 


-3.511 


24.515 


1.267 


24.696 


4.013 


24.531 


0.073 










2000.09.17 


03:27:31 


-3.825 


24.190 


0.742 


29.329 


4.231 


24.343 


1.236 




HS 0002+1635 


23.518 


2.450 


2002.12.02 


01:07:24 


-22.829 


23.518 


2.450 






23.518 


2.450 




WD 0005-163 


15.006 


0.005 


2000.09.16 


03:31:59 


-2.013 


15.057 


1.892 


14.814 


3.648 


15.006 


0.140 










2002.08.04 


10:00:19 


16.515 


15.921 


1.860 


9.051 


4.473 


14.907 


3.445 




WD 0011+000 


25.655 


0.106 


2000.07.14 


07:14:10 


28.209 


23.079 


0.949 


20.950 


2.568 


22.823 


0.978 










2000.07.17 


07:38:21 


27.631 


25.660 


0.657 


25.542 


4.107 


25.657 


0.025 




WD 0013-241 


15.760 


0.061 


2000.09.16 


02:44:05 


-3.848 


15.754 


1.063 


15.797 


2.591 


15.760 


0.020 










2000.09.17 


01:52:24 


-4.237 


13.188 


1.268 


10.630 


2.367 


12.617 


1.505 




WD 0016-258 


44.969 


1.523 


2000.09.16 


03:01:00 


-4.332 


45.801 


1.451 






45.801 


1.451 










2000.09.17 


02:09:57 


-4.713 


44.016 


2.194 


39.586 


6.581 


43.573 


1.879 




WD 0016-220 


10.875 


1.715 


2000.09.16 


05:11:37 


-2.989 


12.101 


0.868 


16.054 


1.894 


12.788 


2.117 










2000.09.17 


03:47:05 


-3.294 


10.506 


0.742 


7.857 


1.757 


10.105 


1.343 


Falc 


WD 0017+061 


-1.247 


3.824 


2002.09.26 


07:34:49 


2.674 


-0.139 


2.876 


-7.848 


7.022 


-1.247 


3.824 


WD 0018-339 


30.744 


0.565 


2002.09.15 


02:14:27 


-6.478 


31.118 


1.042 


29.443 


2.256 


30.823 


0.901 


o 








2002.09.18 


02:33:07 


-7.790 


30.220 


1.117 


21.817 


2.406 


28.729 


4.539 


p 


WD 0024-556 


84.029 


2.130 


2000.08.03 


09:18:35 


-1.148 


84.420 


1.490 


78.216 


5.749 


84.029 


2.130 




WD 0029-181 


42.792 


0.136 


2002.09.26 


08:24:34 


-4.938 


40.163 


1.568 


46.906 


1.747 


43.170 


4.740 


B 








2002.09.27 


06:13:41 


-5.196 


42.249 


1.494 


44.187 


2.475 


42.767 


1.212 


CTQ 


HE 0031-5525 


39.229 


6.206 


2001.12.17 


00:53:54 


-25.943 


37.092 


1.627 






37.092 


1.627 


o 








2002.07.27 


06:08:08 


1.501 


49.230 


1.699 


42.597 


4.058 


48.241 


3.342 


k 


HE 0032-2744 


48.911 


3.650 


2002.09.15 


02:59:55 


-2.771 


51.692 


2.379 






51.692 


2.379 








2002.09.18 


03:06:41 


-4.129 


46.515 


2.208 






46.515 


2.208 


o 


WD 0032-175 


31.342 


1.349 


2002.09.18 


03:21:03 


-0.053 


30.550 


0.995 






30.550 


0.995 


era 








2002.09.25 


06:31:34 


-3.823 


32.492 


1.200 






32.492 


1.200 


ornery, . 


WD 0032-177 


9.875 


3.692 


2002.09.18 


03:29:57 


-0.101 


10.312 


1.452 


6.331 


1.284 


8.078 


2.793 








2002.09.25 


06:01:55 


-3.790 


15.685 


1.756 


9.583 


2.500 


13.668 


4.059 


WD 0033+016 


88.746 


2.451 


2002.09.26 


07:57:13 


2.732 


88.746 


2.451 






88.746 


2.451 


MCT 0033-3440 


47.910 


0.286 


2000.09.16 


04:36:44 


-6.179 


52.870 


1.751 






52.870 


1.751 


?? 








2002.08.15 


09:54:03 


6.272 


47.879 


1.316 


48.012 


2.860 


47.902 


0.071 


t 


HE 0043-0318 


67.616 


1.591 


2002.12.02 


01:30:55 


-26.459 


68.302 


0.967 


70.191 


2.082 


68.637 


1.020 










2003.01.16 


01:44:09 


-29.463 


66.376 


1.124 






66.376 


1.124 


liams 


WD 0047-524 


23.980 


0.368 


2002.07.27 


06:21:40 


3.620 


24.533 


0.795 


22.385 


1.627 


24.119 


1.198 








2002.09.14 


03:12:24 


-10.877 


24.175 


0.716 


19.757 


1.671 


23.488 


2.263 


HS 0047+1903 


27.280 


0.926 


2002.09.27 


05:30:35 


10.719 


27.042 


1.183 


29.084 


3.262 


27.280 


0.926 




WD 0048-544 


22.199 


0.822 


2002.07.27 


06:30:11 


2.495 


22.080 


1.035 


26.192 


2.158 


22.850 


2.268 










2002.09.14 


03:21:12 


-11.578 


23.188 


0.995 


20.319 


0.944 


21.679 


2.026 




WD 0048+202 


35.980 


0.587 


2002.09.19 


04:52:42 


15.032 


34.901 


2.313 


36.145 


1.517 


35.771 


0.807 










2002.09.27 


05:16:50 


11.230 


36.496 


1.197 


39.764 


2.555 


37.085 


1.776 










2002.12.02 


02:09:09 


-18.947 


36.612 


1.150 


31.783 


2.525 


35.783 


2.575 




HE 0049-0940 


27.763 


0.192 


2002.09.26 


08:39:23 


0.509 


27.896 


1.001 


26.624 


2.178 


27.675 


0.682 










2002.09.27 


06:27:15 


0.251 


28.308 


0.783 


26.327 


1.736 


27.973 


1.049 




WD 0050-332 


35.597 


1.801 


2002.07.27 


06:44:31 


14.143 


35.763 


4.027 


36.931 


5.721 


36.150 


0.777 










2002.09.15 


03:26:56 


-3.174 


31.345 


2.874 






31.345 


2.874 










2002.09.25 


06:49:51 


-7.716 


32.107 


2.952 






32.107 


2.952 




WD 0052-147 


56.344 


1.331 


2002.09.26 


08:54:09 


-1.164 


58.289 


1.884 






58.289 


1.884 










2002.09.27 


06:41:44 


-1.422 


56.297 


2.405 


54.873 


3.785 


55.888 


0.912 





TABLE 2 

Apparent Velocity Measurements for Thick Disk DAs 





Adopted 


Date 


Time 


LSR 


Ha 




H/3 


Observation 




^app 


d-Dapp 


Correction 


^app 


dfapp 


^app 


flfapp 


^app 


d^app 


Target 


(km s 1 ) 


(km s 1 ) 


(UT) 


(UT) 


(km s- 1 ) 


(km s 1 ) 


(km s 1 ) 


(km s- 1 ) 


(km s- 1 ) 


(km s i ) 


(km s i ) 


WD 0158-227 


-12.279 


2.795 


2002.09.20 


03:40:17 


-5.823 


-10.305 


1.125 






-10.305 


1.125 








2002.09.27 


07:46:45 


-9.161 


-14.259 


1.127 






-14.259 


1.127 


WD 0204-233 


82.384 


0.304 


2000.07.15 


07:45:18 


8.895 


82.192 


0.944 


80.422 


1.662 


81.760 


1.074 








2000.07.17 


08:47:06 


8.750 


82.364 


1.977 


83.186 


5.477 


82.459 


0.371 


WD 0255-705 


47.623 


0.527 


2000.08.03 


09:55:52 


-35.536 


48.181 


1.547 






48.181 


1.547 








2000.08.05 


08:51:01 


-35.778 


47.374 


1.034 






7.374 


1.034 


WD 0352+052 


-86 636 


2.300 


2002.03.01 


01:09:46 


-35.737 


-87.408 


1.401 






-87.408 


1.401 








2002.09.13 


09:27:33 


19.231 


-84.384 


1.248 


-79.498 


2.216 


-83 208 


2.954 


HE 0409-5154 


23.915 


6.616 


2000.09.15 


09:10:28 


-37.306 


22.243 


2.047 


20.299 


1.634 


21.056 


1.340 








2001.09.13 


09:26:42 


-37.020 


30.762 


1.257 


34.543 


2.406 


31.572 


2.194 


HE 0416-1034 


44.087 


0.191 


2000.12.17 


06:13:24 


-38.231 


43.785 


1.509 


45.786 


3.788 


44.059 


0.973 








2001.01.15 


03:03:25 


-48.043 


40.481 


2.058 


47.436 


1.631 


44.754 


4.787 


HE 0452-3444 


-10.387 


0.086 


2000.12.13 


06:27:39 


-44.215 


-9.733 


1.464 


-13.936 


3.554 


-10.343 


2.094 








2001.01.15 


02:44:23 


-51.955 


-9.238 


2.210 


-13.956 


3.712 


-10.472 


2.932 


HE 0508-2343 


79.982 


2.028 


2001.04.07 


00:23:02 


-59.490 


79.124 


1.786 


80.903 


1.307 


80 283 


1.198 








2001.04.09 


00:48:48 


-59.118 


76.884 


1.817 


68.789 


1.966 


73.155 


5.706 


WD 0732-427 


36 295 


0.558 


2001.04.09 


01:16:13 


-68.414 


36 003 


0.756 


37.866 


1.433 


36.409 


1.087 








2001.05.03 


23:55:10 


-69.286 


33.730 


1.190 


40.950 


2.679 


34.920 


3.788 


HS 0820+2503 


41.390 


4.704 


2003.02.18 


02:42:33 


-33.578 


41.390 


4.704 






41.390 


4.704 


HE 1124+0144 


46.460 


0.144 


2000.07.01 


23:28:45 


-55.699 


45.879 


1.394 


48.423 


2.799 


46.384 


1.435 








2000.07.02 


23:34:25 


-55.550 


47.083 


1.185 


42.748 


3.337 


46.598 


1.932 


WD 1152-287 


51.828 


1.745 


2000.07.11 


00:57:00 


-70.466 


50.760 


1.594 






50.760 


1.594 








2000.07.14 


23:01:27 


-70.053 


53.254 


1.841 






53.254 


1.841 


WD 1323-514 


-40.353 


0.770 


2001.05.15 


01:29:38 


-42.147 


-40.093 


0.793 


-41.065 


2.400 


-40.189 


0.409 








2001.06.08 


02:12:18 


-50.778 


-42.707 


1.018 


-40.445 


1.821 


-42.168 


1.362 


WD 1334-678 


22.946 


0.649 


2000.07.30 


00:47:19 


-58.211 


22.384 


1.473 


24.018 


3.072 


22 689 


0.900 








2001.05.15 


01:42:37 


-40.116 


24.128 


1.309 


20.157 


4.143 


23.767 


1.612 


HS 1338+0807 


64.191 


2.405 


2001.08.18 


23:52:07 


-28.395 


64.191 


2.405 






64.191 


2.405 


WD 1410+168 


15.184 


2.748 


2002.04.23 


06:08:33 


7.586 


15.753 


0.911 


8.557 


3.107 


15.184 


2.748 


WD 1426-276 


55 800 


0.559 


2000.07.05 


03:56:09 


-43.083 


55.742 


1.182 


55.448 


3.136 


55.705 


0.137 








2000.07.06 


02:38:26 


-43.227 


57.693 


1.045 


56.772 


1.755 


57.452 


0.572 


HS 1432+1441 


84.966 


3.366 


2001.08.16 


00:01:42 


-14.425 


82 686 


1.642 


77.964 


3.260 


81.730 


2.683 








2001.08.21 


00:15:53 


-13.555 


85 505 


1.995 


88.323 


2.297 


86.717 


1.973 


WD 1507+021 


50.907 


3.076 


2002.06.18 


01:25:19 


-7.271 


52.338 


2.198 


47.601 


3.340 


50.907 


3.076 


WD 1614-128 


95 682 


0.052 


2000.06.06 


06:16:38 


5.840 


94.870 


1.239 


99.368 


2.581 


95.713 


2.482 








2000.06.08 


02:20:56 


5.346 


94.931 


1.088 


102.053 


3.282 


95.637 


3.009 


WD 1716+020 


10.643 


1.291 


2002.04.23 


09:03:16 


48.550 


11.038 


0.846 


8.531 


1.956 


10.643 


1.291 


WD 1834-781 


60.781 


0.111 


2000.07.06 


04:36:43 


-35.096 


60.857 


0.919 


60.355 


2.395 


60.792 


0.237 








2000.07.13 


05:04:33 


-37.059 


59.611 


1.028 


62.217 


1.793 


60.255 


1.589 


WD 1952-206 


62.701 


0.113 


2000.07.06 


05:09:15 


29.248 


62.167 


0.870 


63.309 


2.141 


62.329 


0.563 








2000.07.13 


04:43:35 


25.831 


62.683 


0.956 


62.931 


2.356 


62.718 


0.122 


WD 2029+183 


-97.198 


1.373 


2002.04.24 


09:34:58 


74.163 


-98.949 


2.029 


-96.770 


1.646 


-97.635 


1.507 








2002.08.05 


04:50:45 


51.952 


-93.455 


1.652 


-98.577 


2.467 


-95.041 


3.349 


WD 2322-181 


39.763 


1.328 


2000.07.13 


06:16:26 


38.666 


42.172 


1.462 


38.142 


2.548 


41.173 


2.459 








2000.07.16 


05:46:18 


37.813 


39.731 


1.648 


36.877 


3.217 


39.138 


1.638 


WD 2350-083 


84.801 


0.992 


2002.07.11 


09:56:33 


49.133 


84.652 


1.390 


85.845 


2.341 


84.963 


0.740 








2002.09.13 


07:51:37 


24.339 


80.737 


1.489 


86.827 


3.323 


81.756 


3.214 
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TABLE 3 

Mean Apparent Velocities 



Sample 


# of WDs 


(l>app> 

(km s ) 


<5(Uapp) 

(km s ) 


(km s ) 


(km s — 1 ) 


(M/R) 
(M /R ) 


5{M/R) 
(M @ /R©) 


Main 


449 


32.57 


1.17 


24.84 


1.51 


51.19 


1.84 


Ho 


449 


32.69 


1.18 


24.87 


1.78 


51.37 


1.85 


H/9 


372 


31.47 


1.32 


25.52 


3.17 


49.45 


2.07 


Thick 


26 


32.90 


9.59 


48.99 


1.57 


51.70 


15.07 



TABLE 4 

Mean Masses 



Sample 


# of WDs 


(km s — ) 


S{va,pp) 
(km s — 1 ) 


<Teff> 

(K) 


(K) 


(M Q ) 


Main 


449 


32.57 


1.17 


19400 


9950 


0.647 


Thick 


26 


32.90 


9.59 


19960 


11060 


0.652 


Hot a 


366 


31.61 


1.22 


21670 


9700 


0.640 


Cool a 


75 


37.50 


3.59 


9950 


1090 


0.686 



S(M) 
(M s ) 



+0.013 
-0.014 
+0.097 
-0.119 
0.014 
+0.035 
-0.039 



"Hot" refers to WDs with T cft > 12000 K and "cool" to WDs with 12000 K > T cff > 7000 K. 



TABLE 5 

Mean DA Masses From Selected Previous Studies 





(M) 


5 (A*) 


cm 


Assumed 




Study 


# of WDs (M Q ) 


(M ) 


(M ) 


Method H-Layer a 


Notes 



Kocstcr ct al. (1979) 
Koester (1987) 
McMahan" 7l989) 
Weener fc Reid (1991) 
Bergeron et al. (1992) 
Bragaglia ct al. (1995) 
Bergeron et al. (1995b) 

Reid (1996) 
Vennes et al. (1997) 
Finlev et al. (1997T 

Silvestri et al. (2001) 
Madei et al~T2004) 
Liebert et al. (2005) 

Kepler et al. (2007) 
Tremblav JTBergeron (2009) 

Koester et al. (2009b) c 
~ (2009b) c 

(2009b) overlap 



Koester et a 



Koester et 



This Work 



122 


0.58 


0.10 





12 b 


Photo 


Thin/No 




9 


0.58 







11 


GRS 


Thin/No 


CPM WDs 


50 


0.523 


0.014 






Spectro 


Thin/No 




35 


0.63 


0.03 






GRS 


Thin/No 


CPM WDs 


129 


0.562 







137 


Spectro 


Thin/No 


T aH > 14000 K 


42 


0.609 







157 


Spectro 


Thin/No 


T eff > 12000 K 


129 


0.590 







134 


Spectro 


Thick 


Revised Bergeron et al. (19921 
w/ thick H-layers 


34 


0.583 


0.078 






GRS 


Thick 


CPM WDs 


110 


0.56* 








Spectro 


Thin/No 


75000 K > T cS > 25000 K 


174 


0.570* 







060* 


Spectro 


Thick 


T a « > 25000 K 
some w/ cool companions 


41 


0.68 


0.04 






GRS 


Thick 


CPM WDs 


1175 


0.562* 








Spectro 


Thick 


T e g > 12000 K 


298 


0.603 
0.572* 








134 
188 


Spectro 


Thick 


T cff > 13000 K 


1859 


0.593 


0.016 






Spectro 


Thick 


T cS > 12000 K 


- 250 


0.649 








Spectro 


Thick 


40000 K > T cfi > 12000 K 
overlap w/ Liebert et al. 


606 d 


0.567 c 


0.002 e 





142 c 


Spectro 


Thick 


SPY 


512 d 


0.580 c 


0.002 c 





136 e 






SPY non-binary thin disk WDs 


441 


0.575 e 


0.002° 





128 e 








449 


0.647 


+0.013 
-0.014 






GRS 


Thick 


SPY non-binary thin disk WDs 



Note. — Masses marked with an asterisk arc peaks/widths of mass distributions from Gaussian fitting. 
a Hydrogen layer mass used in mass-radius relation from evolutionary models. "Thick" corresponds to M^/M* ~ 10 -4 and "Thin/No" 
A / 1 1 / M k < 10 — 8 or no hydrogen layer, 
k Two-thirds of the stars arc within 0.12 Mq. 

C Masses do not appear in this reference. We compute masses from the published values of log g and T e ff using the mass-radius relation fr 
evolutionary models. 

Excludes double degenerates. 
e We compute these means, uncertainties, and standard deviations (sec Scction l5.2|> . 



